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MEASUREMENTS OF AVERAGE HEAT-TRANSFER AND FRICTION COEFFICIENTS 
FOR SUBSONIC FLOW OF AIR IN SMOOTH TUBES AT HIGH 
SURFACE AND FLUID TEMPERATURES ^ 

By Leeot V. Humble, Wahren H. Lowdehmilk, and Lelaks C. Desmon" 


SUMMARY 

An investigation of forced-convection, heat transfer and 
associated fressvxe drops was conducted with air flowing 
through smooth tubes for an over-all range of surface tempera- 
ture from BS5° to 3050° R, inlet-air temperature from 535° 
to 1500° R, Reynolds number up to 500,000, exit Mach number 
up to 1, heat flux up to 150,000 Btu per hour per square foot, 
length-diameter ratio from 30 to 120, and three entrance con- 
figurations. Most of the data are for heat addition to the air; 
a few results are included for cooling of the air. The over-all 
range of surface-to-air temperature ratio was from 0J[6 to 3.5. 

Correlation of the measured average heat-tran^er and friction 
coefficients with heat addition by conventional methods wherein 
the physical properties of the air were evaluated at the average 
air temperature resulted in considerable decreases in both die 
NusseU number and the friction coefficient at constant Reynolds 
number in the turbulent region, as the ratio of surface tempera- 
ture to air temperature was increased. The effect of sufface- 
to-air temperature ratio was eliminated by evaluating the 
physical properties of air, including density, in the Reynolds 
number, the Nusselt number, and the friction coefficient at a 
temperature, higher than the average air temperature. 

Correlation of the heat-transfer coefficients was also affected 
by increases in the average air temperature, produced by in- 
creasing the inlet-air temperature, which resulted in a decrease 
in the Nusselt number. This effect of temperature level was 
reduced by using the smallest variation of thermal conductivity 
with temperature available in the literature and was eliminated 
when- the thermal conductivity was arbitrarily assumed to 
■vary as the square root of temperature. 

INTRODUCTION 

A large amount of data is available in tlie literature on 
forced-convection heat transfer from sinfaces to fluids. 
Most of these data, however, have been obtained at relatively 
low surface temperatures and heat-flux densities, and do not 
extend into the range of high temperature and flux that is 
of interest in many current engineering applications. In- 
asmuch as convective heat transfer is a boundary-layer 
phenomenon, the severe temperature and velocily gradients 


in the fluid film adjacent to the surface (and attendant 
variation in fluid properties), which accompany heat transfer 
at high flux densities, make extrapolation of existing data 
for low flux densities uncertain. This possibility has been 
previously recognized; for example, in references 1 and 2 
it is indicated that heat transfer and friction may depend 
on both surface and fluid temperatures. 

An experimental investigation was undertaken at the 
NACA Lewis laboratory during 1948-50 to obtain heat- 
transfer and related pressure-drop information for a wide 
range of surface and fluid temperatures and heat flux. As 
part of the general program, an investigation was made 
with air flowing through smooth tubes. The effects of such 
variables as surface temperature, inlet-air temperature, and 
tube-entrance configuration on heat transfer and pressure 
drop have been investigated and are reported in references 
3 to 6. The results of references 3 to 6 are summarized 
herein and, in addition, previously unpublished data showing 
the effect of tube length-diameter ratio and inlet-air tempera- 
ture on heat-transfer and friction coefficients at high surface 
temperature and heat flux are presented. Most of the data 
are for heat addition to the air; however, a few previously 
unpublished results are included for heat extraction from 
the air. 

APPARATUS 


AKEANGEMENT 


A schematic diagram of the equipment used in the investi- 
gation is shown in figure 1 (a) . Compressed air was supplied 
through a pressure-regulating valve, cleaner, and surge tank 
to a second pressme-r^ulating valve where the flow rate was 
controlled. From this second valve, the air flowed through 
metering devices and a preheater into a mixiTig tank, which 
consisted of three concentric passages so arranged that the 
air made three passes through the tank before entering the 
test section. Baflles were provided in the central passage 
to insure thorough mixing of the air before it entered the 
test section. From the test section the air flowed through a 
second mixing tank and was discharged to the atmosphere. 
The tffit section, mixTug tanks, and adjoining piping were 
thermally insulated. 


1 Supersedes NAQA. RM E7L31, “Heat Transfer Itozn HJgh-Temperatnre Surfaces to PlnMs. I — Preliminary fiivestfeation with Air fn Inconel Tube ^wftli Rounded Entrance > 
Inside Diameter of0.4lncli, and Length of 24 Riches" by Leroy V. Humble, ‘Warren H. Lowdsnnfik, and Jdilton Grele, 1948; XACA RM E8L08, “Heat Transfer frwn Hi^Temperature 
Sorfaces to Flnlds. H— Correlation of Heat-Transfer and Friction Data for Air Flowing In Inconel Tnbe with Rounded Entrance" by Warmi H. Lo^dermHk and Mfltoai D. Grele, 1949; 
NACA RM E80E23, “Rifluence of Tnbe-Entrance Conflguratton on Average Heat-Transfer Coefficients and Friction Factocs for Air Flowing In an Inconel Tube" by Warren H.LowdermiQ: 
and MUlon D. Grele, 1950; NACA RM E60H23, ‘^Correlation of FOTced-Couvectlon Heat-Transfer Data for Air Flowing in Smooth Platinnm Tube with Long-Approach Entrance 
at High Surface and Inlet-Alr Temperatures" by Leland G. Beamon and Eldon W. Rattis, i 65 q. 
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Votfrmafer 



(a) Schematic diagram of setup. 



(b) Tjiplcal test section showing thermocouple and pressure-tap locations. 
FiaUKi 1.— Arrangement of apparatus Ibr heating air. 


The temperature of the air entering and leaving the test 
section was measured by thermocouples located downstream 
of the mixing baffles in the entrance and exit mixing tanks, 
respectively. 

For the runs ivith heat addition to the air, electric power 
was supplied to the test section from a 208-volt, 60-cycle 
supply line through an autotransformer and a step-down 
power transformer. The low-voltage leads from the power 
transformer were connected to flang;es on the test section by 
flexible copper cables. A voltmeter, an ammeter, and a 
wattmeter were provided to measure the electric-power 
input to the test section. 

For the runs with heat extraction from the air, the electric- 
power supply was disconnected and the test section was sur- 
i-ounded with a water jacket. The warmed air from the pre- 
heater was cooled in the water-jacketed test seption. 

TEST SECTIONS 

Most of the investigation was conducted using test sections 
fabricated from commercial Inconel tubing. The range of 
tube-wall temperature was extended beyond that attainable 
with Inconel by. using one teat section made of platinum 
(reference 6), which wa^ installed in a separate setup similar 
to that shown in figure 1 (a). One of the test sections used 





(а) Long approach. 

(б) Bigbt-angle e<]ge. 

( 0 ) Bcllmoutb. 

Fig CEE 2.— EntranccB to test sections . 

for heating the air is illustrated schematically in figure 1 (b). 
Flanges were attached to both ends of the tube to provide 
electric contact with the power supply. The dimensions of 
the various tubes used as test sections in the investigation 
were as follows; 


Tube material 

Inside 

diameter 

(In.) 

Outside 

diameter 

(in.) 

Heat- 
transfer 
length 
. (In.) 

Length- 

diameter 

ratio 

Incone] — 

0.402 

0,500 

12 

30 

Do 

.402 

.wo 

24 

60 

Do 

,402 

.500 

48 

120 

Platinum - 

. .£2f 

.OSS 

24 

40 


Outside-wall temperatures wei’e measured at a number of 
stations along the tube by means of thermocouples and self- 
balancing indicating-type potentiometere. Chromel-alumel 
thermocouples were used with the Inconel test sections, and 
platinum — platinum-rhodium thermocouples, tvith the plati- 
num test section. In general, at each station along the tube, 
two thermocouples were located diametrically opposite and 
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the average of the two temperatures was taken as the outside- 
wall temperature at the station. Static-pressure taps were 
located at mtervals along each of the Inconel test sections. 

Investigations were made with heat addition to the air 
using three types of entrance on the 24-inch Inconel tube 
(reference 5). The entrances, shown in figure 2, were 
(a) long approach, (b) right-angle edge, and (c) bellmouth. 

The cooling data were obtained with an Inconel test section 
having a length-diameter ratio of 60 and a bellmouth en- 
trance. 


PROCEDURE 


The following general procedure was used in obtaining the 
experimental data: For the heating runs, the inlet-air tem- 
perature was set at the desired value, the inlel-air pressure 
was adjusted to give the rnTnirnmn desired flow rate (mini- 
mum Reynolds number), and the electric-power input was 
adjusted to give the desired tube-waU. temperature. After 
equilibrium conditions had been attained, electric-power 
input, flow rate, temperatures, and pressures were recorded. 
The air flow and the power input were then increased in 
increments (to increase the Reynolds number while main- 
taining constant wall temperature) and data were recorded 
after each incremental increase. The foregoing procedure 
was repeated for a range of tube-waU temperature with each 
of the various test sections and entrance configurations. 

For the cooling runs, the Met-air pressure was adjusted to 
give the minimum desired air-flowrate at the desired inlet-air 
temperature. The water-flow rate in the water Jacket was 
set at a constant value. After equilibrium conditions w'ere 
obtained, the data were recorded aS in the heatmg runs. 
The inlet-air pressm’e was then incrementally varied to cover 
a range of Reynolds number while maintaining constant 
inlet-air temperature. The foregoing procedure was repeated 
for a range of inlet-air temperature. 

The over-all range of conditions for which data were ob- 
tained is summarized in the following table : 


. 

Refer- 

ence 

Tube ma- 
terial 

Length- 

diameter 

ratio 

Entrance 

type 

Maxi- 

mom. 

bollr' 

Reynolds 

nnmber 

Average 

tuhe-wall 

tempera- 

ture 

CH) 

■■ 

Inlet-air 

tempera- 

ture 

UR) 

MaxL 

mnm 

exit 
A'Ia(^ 
nmnfaer | 

3 

Inconel 

60 

Bellmouth-- 

260,000 

6S0-1700 

535 

LO 

4 


GO 


500.000 

375.000 

535-2050 

535 

LO 

6 

-_.do 

60 

Long ap- 
proacot 

rffrht ftnelfi 

635-1960 

535 

I.O 

r G 

Platlmmi- 

46 

Long ap- 
proach 
Bellmouth. . 

320,000 

980-3050 

535-1160 

0.6 


Inconel 

30 

800,000 

635-1685 

635 

1.0 


120 


282,000 

250,000 

635-1850 

535-1460 

LO 


. do 

60 

do„ 

585-585 

585-1600 

LO 


« Previously unpablfdied, 
b Cooling data. 


The over-all range of heat-flux density encountered in the 
investigation was from 500 to 150,000 Btu per hour per 
square foot of heat-transfer area. 

21363T— 58 ^28 


SYMBOLS 


The following symbols are used in the report: 

Cj specific heat of air at constant pressure 

(Btu/ab) (°F)) 

D insi de diameter of test section (ft) 

f average friction coefl5.cient 

j(^, f, modified average friction coefficient 

O mass velocity (mass flow per unit cross- 

sectional area) (lb/(hr) (sq ft)) 
g acceleration due to gravity (4.17X10® 


A 

k 


kc 


L 

P 

P 

Ap 


Apfr 

Q 

R 

r 

S 

T 

T, 

Tf 

Tojs 

To 

T. 

t 

V 
Va, 
W 

X 

Y 
/t 
P 

Pap 


ft/hr") 

average heat-transfer coefficient (Btu/(hr) 
(sq ft) (°F)) 

thermal conductivity of air (Btu/(hr) 
(sq ft) (°F/ft)) 

thermal conductivity of test-section ma- 
terial (Btu/(hr) (sq ft) (°F/ft)) 
heat-transfer length of test section (ft) 
absolute -total or stagnation pressure 
ab/(sq ft)) 

absolute static pressure (lb/(sq ft)) 
over-all static-pressure drop across test 
section (lb/(sq ft)) 

friction static-pressure drop across test 
section (lb/(sq ft)) 
rate of heat transfer to air (Btu/hr) 
gas constant for air (53.35 ft-lb/(lb) (“F)) 
radius of test section (ft) 
heat-transfer area of test section (sq ft) 
total or stagnation temperature (°R) 
average bulk temperature, defined by 
(Ti-FT*)/2, C>R) 

average film temperature, defined by 
(r,-hT,)/2, (°R) 

average film temperature, defined by 
T,+0.75 (T.-T,), (°R) 
average outside-waU. temperature of test 
section (°R) 

average inside-wall (surface) temperature 
of test section (°R) 

.static temperature (°R) 
velocity (ft/hr) 

average velocity, defined by G/Pav> (ft/hr) 
air flow (Ib/hr) _ 

distance from entrance of test section (ft) 
ratio of specific heats of air 
absolute viscosity of air (lb/(hr) (ft)) 
density of air (lb/(cu ft)) 
average density of air defined by (pi+Pif) 
R ik+fi), (lb/(cu ft)) 


Cpfifk 

ODIf. 

PfViDffif 

Pojsl' l)DIP0.7t 

hDjk 


Prandtl number 
! Ie 3 Tiolds number 

modified Reynolds number 

Nusselt number 
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test-section entrance 
test-section exit 
inner surface of test section 
outer surface of test section 
bulk (when applied to properties, in- 
dicates evaluation at average bulk 
temperature T^) 

film (when applied to properties, in- 
dicates evaluation at average film 
temperature TV) 

film (when applied to properties, in- 
dicates evaluation at 
surface (when applied to properties, 
indicates evaluation at average inside- 
wall temperature of test section T,) 

METHOD OF CALCULATION 

Physical properties of air. — The physical properties of ab- 
used in the present investigation are shown in figure 3 as 
functions of temperature. The solid curves are from references 
7 and 8, and are in agreement with the values of reference 9. 
The dashed curve for thermal conductivity is from reference 
10. The dotted curve represents values of thermal con- 
ductivity assumed herein to vary as the square root of tem- 
perature and its use will be discussed in a subsequent section. 
The thermal conductivities from references 7 to 10 were 



Temperafura, T, “R 

Ftiiuni 3.— Variation of thermal oonduotlvlty, absolute vlaoislty, apeolhc heat, and Prandtl 
number of air with temperature. 


Subscripts: 

1 

2 

i 

0 

b 

/ 

0.75 

< 


extrapolated considerably beyond the range of expi'rimental 
data in the reference reports and have been extrajiolatcd 
still further in figure 3 to cover the temperature range of the 
present investigation. The upper temperature limits of the 
experimental data upon which values of thermal conductivity 
from references 7, 8, and 10 are based are 1050°, 1300°, and 
1215° R, respectively. 

Except where other\vise stated, the physical property data 
of references 7 and 8 arc used in the present report. 

Determination of surface temperature, — The average 
outside-wall temperature of the test section T, was obtained 
by measuring the area under a curve of the tempei-atiire 
distribution (as obtained from thermocouple readings) aicng 
the heat-transfer length and dividing the area by this length. 

For the heating runs, the average inside-wall (surface) 
temperature T, was then calculated by the following equation 
(derived in reference 11); 





( 1 ) 


In equation (1), the assumptions are made that heat is gen- 
erated uniformly across the tube-wall thiclvncss and that the 
heat flow from every point in the tube material is radially 
inward. For the present investigation, the radial temperature 
drop tlirough the tube wall, calculated from equation (1), 
was very small compared with the difference between the 
average inside-wall and the air bulk temperature. 

For the cooling runs, the thermocouples were imbedded in 
the test-section wall and were assumed to measure the inside- 
wall temperature, inasmuch as the calculated temperature 
drops tlu-ough the wall were small. 

Heat-transfer coefficients. — The average heat-transfer 
coeffici^t h was computed from the expei-imcntal data by 
tlie relation * 


Trc,.,(T,-r.) 


( 2 ) 


The use of the total air temperature Tt, instead of adiabatic 
wall temperature was considered justified because the differ- 
ence between these temperatures was small compared with 
tlie difference between the average surface and air (or 
adiabatic wall) temperatures. 

Correlation of the average heat-transfer coefficients with 
the pertinent variables is discussed in the section RESULTS 
AND DISCUSSION. 

Friction coefficients. — Friction data were obtained witii 
and without heat transfer. Average friction coefficients were 
calculated from the experimental pressure-drop data as fol- 
lows: ff'he friction pressure drop was obtained by subtracting 
the calculated momentum pressure drop fi-om the measui-od 
static-pressure drop across the test section. Thus 
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where U and t 2 , the absolute static temperatures at entrance 
and exit of the test section, respectiTelj, were, in general, 
calcidated from the measured values of air flow, static pres- 
sm-e, and total temperature by the following equation, which 
is obtained by combining the perfect gas law, the equation of 
continuity, and the enei^ equation: 


yff (P\- 1 tv Tg fpYJ . oy - yg /yV 
'(y-l)R\Gj ■^YL(r-l)BW J ^ (7-l)SW 


(4) 


For the bellmouth entrance, it was found that the entrance 
static temperature U could be represented with sufficient 
accuracy by the relation 


7-1 



In equation (5) the total pressure at the test-section 
entrance Pi was assumed to be equal to the static pressure in 
the entrance mixing tank, where the velocity was negligible. 

Average friction coefidcients were calculated by dividing 
the friction pressure drop by four times the length-diameter 
ratio of the tube times a dynamic pressure. The density in 
the dynamic pressure was computed in three ways, resulting 
in three deflnitions of friction coefficient. 

The first of these friction coefficients was based on a 
density evaluated at the average static pressure and tem- 
perature of the air in the tube; thus 


where 


jf ^'Pfr 9 Pai^Pfr 

j A Pat Vat^ o -Li 
2g 

(6) 

„ _ I f Pi+P^\ 

** 

(7) 


In the other two definitions of friction coefficient, the 
density was evaluated at^the film temperature Tf and at 
the average surface temperature T„ resulting in equations 
(8) and (9), respectively: 


//=- 


/.=- 


^Pfr fPoA r /2TA 

L PrVJ \Pf)^ Kti^tJJ 
D 2g 


( 8 ) 


as \Pt/ 


D 2g 

Correlation of the friction data is discussed in the following 
section. 

RESULTS AND DISCUSSION 

AXIAL WALL-TEMPERATUBE DISTEmUTIONS 


Representative axial waU-temperature distributions 
obtained with heat addition to the air in an Inconel tube 
having a length-diameter ratio of 60 and a bellmouth entrance 


are shown in figure 4. The outside-waU temperature is 
plotted against the ratio of the distance from the tube 
entrance to the heat-transfer length x/L for five different 
rates of heat transfer to the air. The rate of heat transfer 
to the air, the air flow, the average inside tube-wall tem- 
perature, and the temperature rise of the air are tabulated. 

The wall temperature increases almost linearly over the 
major portion of the tube length, and the slope of this part 
of the curve increases with an increase in rate of heat input. 
The large axial-temperature gradients at the entrance and 
exit of the test section are the result of conduction losses 
through the electric connector flanges and cables. 

The curves in figure 4 are typical; s imil ar trends were 
obtained with other entrances and for other length-diameter 
ratios. 



t.engfh of fesf section. ^ ’ L 

Figure 4.— Representative outsfde'wall temperature dlstrfbiitioii for various nmniintg of heat 
input to air. BeTnolds nnmber, approximatel? 140,000; Inconel tube; lengtbrfUameter 
ratio, 60; bellmoutli entrance; inlet-air temperature, 58S° R. 


COBBELATION OF HBAT-TBANSFEB COEFFICIENTS 


Conventional correlation based on. average air tempera- 
ture. — ^Forced-convection, turbulent-flow, heat-transfer co- 
efficients are generally correlated with other pertinent vari- 
ables by means of the familiar relation in which Nusselt 
number divided by Prandtl number to a power (generally 
hD l/c 

0.4) ) is plotted against Eeynolife number GD/fi, 


and in which the physical properties of the fluid (specific 
heat, viscosity, and thermal conductivity) are usually eval- 
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uatcd at the average fluid temperature. Results obtained 
by beating air flowing tbrougb an Inconel tube having a 
length-diameter ratio of 60 and a belhnouth entrance 
(reference. 4) are plotted in this manner in figure 5 (a). 
The data shown are for an inlet-air total temperature of 
535° R and cover a range of Reynolds number from about 
7000 to 600,000. Included for comparison is the average 
line obtained by McAdams (reference 2) from correlation 
of the results of various investigators. The equation cor- 
responding to this reference line is 



(a) Physical properties of air evaluated at averaee bulk temperature. 

(b) Physical properties of air evaluated at average surlhoe temperature. 

Figum J. — Conventional methods of correlatinst heatdransfer ooeffldents. Inconel tuba; 

lengttwllameter ratio, 60; bellmouth entrance; Inlet-alr temperature, 536* R. 


A family of parallel lines representing the various siirfaec- 
to-air temperature ratios and having slopes of about 0.8 at 
Reynolds numbers above approximately 25,000 is obtained. 

The data for low surface temperatures (T„ 605® R) arc in 
reasonable agreement with the reference line; however, as 
the surface temperature is raised the tlata fall progressively 
below tlie reference line. For example, at a Reynolds 
number of 70,000, an increase in T,/Ti, from 1.1 to 2.5 
(corresponding to an increase in surface temperature from 
605° to 2060° R) results in a decrease in the onlinate 
parameter of about 38 percent. Inasmuch as the corre- 
sponding change in Prandtl number to the 0.4 power is 
small, about the same percentage decreaso occurs in Nusselt 
number. 

The results shown in figure 5 (a) are typical of those 
obtained with other entrance configurations and length- 
diameter ratios. 

Conventiosal correlation based on surface temperature. - 
The heat-transfer data of figure 5 (a) are replotted in fig- 
ure 5(b), wherein the physical pi-opertics of the air are evaluated 
at the average inside-wall temperature T,. The line repre- 
senting equation (10),' but W'ith the physical properties 
evaluated at T„ is included for (fomparison. The data again 
show a separation with surface temperature level, which in 
this case is greater tlian when the proptTties are evaluated 
at the .bulk temperature. The use of any temperature 
between these extremes would, of course, cause an inter- 
mediate amount of separation of the data. 

In references 1 and 12, heat-transfer coefficients were 
predicted for fully developed flow of fluids having a Pi’andtl 
number of 1.0 in smooth tubes. The analyses took into 
consideration the radial variation of fluid pi'operties and 
indicated an effect of surface-to-fluid temperature ratio 
analogous to that shown in figure 5. 

Effect of Burface-to-fluid temperature ratio.— The magni- 
tude of the effect on the heat-transfer correlations of surface- 
to-fluid temperature ratios gi’eater than 1.0 is 'illusti’atcd 
more clearly in figure 6. Liuc.s OiTo shown for which the air 
properties are evaluated at bulk temperature T*, film temper- 
atm-e Tf, and surface temperature T,. The lines correspond- 
ing to evaluation of properties at Tt, ami T, are cross plots 
of the data at a Reynolds number of 100,000 from fig- 
ures 5 (a) and 5 (b), respectively. The lino corresponding to 
evaluation of air properties at Tf was obtained from a plot 
similar to those of figure 5. 

A series of straight lines is obtained, W'hich intei'scct at 
T,ITi of 1.0, at which point tlie ordinate parameter has the 
value predicted by equation (10). The lines have increas- 
ingly greater negative slopes as the temperature at which 
the fluid properties are evaluated is increased. For e.xample, 
an increase in T./Tj from 1.0 to 2.5 results in a decrease in 
the ordinate parameter of about 40 percent when the prop- 
erties are evaluated at Tt, of 48 percent when evaluated at 
Tf, and of 64 percent when evaluated at T,. 
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Appreciable en’or can be made (fig. 6) if conventional heat- 
transfer correlation equations' are used to compute heat- 
transfer coeJEcients for surface-to-fluid temperature ratios 
appreciably greater than l.O. 

Modified correlation based on surface temperature. — 
The slope of the line in figure 6 corresjwnding to evaluation 
of the fluid properties at T, is —0.8. The ordinate in this 
case is {kDlk,)l(cpjii,fkf)° and the corresponding Keynolds 
number is defined by GDffi,. Inasmuch as at T,fTt=l.Q 
equation (10) applies, it might be expected that the effect 
of T,/Tj could be eliminated by incorporating the term 
(T’s/T,)®’® in equation (10) with the properties evaluated at 
T,. Thus 




( 11 ) 


number divided by Prandtl number to the OA power is 
plotted against the modified Reynolds number 

The use of the modified Reynolds number eliminates the 
trends with respect to surface-to-bulk temperature ratio, so 
that for the range of surface temperature and Reynolds 
number shown the data can be represented with good accu- 
racy by a single line. For modified Reynolds numbers 
above about 10,000, the data can be represented very well 
by equation (14). 

Effect of entrance configuration. — ^Heating data obtained 
with an Inconel tube having a length-diameter ratio of 60 
and both long-approach and rig^t -angle-edge entrances are 
shown in figure 8. The coordinates are the same as in fig- 
ure 7, and the line representing the data of figure 7 (beUmouth 
entrance) is included in the figure. 


If the difference between total and static temperatures is 
neglected. 


T, pt 


( 12 ) 


and the following relation is obtained by letting 6=PbV„i 


(f)(0=^ 


(13) 


Equation (13) defines a modified form of Reynolds number 
in which the conventional mass velocity Q is replaced by the 
product of air density evaluated at the average surface tem- 
perature p, and air velocity evaluated at the average bulk 
temperature T'j. Substituting equation (13) in equation (11) 
gives 

f =0.0.3(^X'(^r (14) 

The data of figure 5 are replotted in figure 7. Nusselt 



FiGCTE 6.— Variation, of Nusselt number dlTided by Prandtl munber to 0.4 power wIHi ratio 
of soi&oe-to-bnlk temperature uslng^ conrentlanal methods of correlathiK heat-transler 
coefficients. Reynolds number, 100,000; Ihoond tube; lengtltdiameter ratio, 00; bellmoath 
entranee; Inlet-air temperature, 030° E. 



Modified Reynolds number, 


PiouBE 7. — Modified method of correlating beat-transfer ooefficlents. Inconel tnbe; length- 
diameter ratio, 60; beUmouth entrance; Inlet-efr temperature, 035° B; physical properties of 
air eralnated at snriace temperature. 


The long-approach and beUmouth entrances tend to pro- 
vide fuUy developed and imiform velocity profiles, respec- 
tively, at the inlet to tike test section. The right-angle-edge 
entrance is typical of most tube-type heat exchangers. 

As shown in figure 8, the data obtained with the long- 
approach and right-angle-edge entrances can be represented 
with good accuracy by the line which best fitted the data for 
the beUmouth entrance shown in figure 7. It may be con- 
cluded that the average heat-transfer coeflicient is negligibly 
affected by the entrance configurations investigated, at least 
for length-diameter ratios of the order of 60 or greater. 

Effect of length-diameter ratio. — ^Heating data obtained 
with Inconel tubes having length-diameter ratios of 30 and 
120 and with a platinum tube having a length-diameter 
ratio of 46 (reference 6) are shown in figure 9 (a). The 
coordinates are the same as in figures 7 and 8 and the line 
(shown solid) representing the data of figures 7 and 8, which 
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are for a length-diameter ratio of 60, is included. The In- 
conel tubes had bellmouth entrances, and the platinum tube 
had a long-approaeh entrance ; however, as indicated in figure 
8, the effect of entrance configuration is negligible. 

A small trend with respect to Icngth-diametei' ratio is 
evident in tJie turbulent region. The data for length- 
diameter ratios less than 60 tend to fall parallel to, but above, 
the reference line, whereas the data for a length-diameter 
ratio of 120 fall below the line. This effect of length- 



Figube S.^Effect of entranoo ocmfiguratloa tialng modified method of correlatlog heat* 
transfer ooeffidents. Inoonel tube; length-diameter ratio, €0; Inlet-air temperatore, B; 
phyilcal properties of 9\x evaluated at siirfocie temperature. 


diameter ratio has been observed by other investigators (see, 
for example, references 1 and 13) and was accounted for by 
including a power function of the length-diameter ratio in 
the general correlation equation. The exponent 0.054 is 
quoted in reference 1 and tlie value 0.1 is quoted in refer- 
ence 13. The exponent 0.1 applies best for the present data 
and accordingly the data of figures 7 and 9 (a), corrected for 
the effect of length-diameter ratio, arc replotted in figure 9 (b). 

The correction for length-diameter ratio appreciably im- 
proves the correlation for modified Reynolds numbers above 
approximately 10,000, and in this region the data can be 
represented by the following equation: 


kD 


0.034 


^rmxir 


(15) 


The constant 0.034 is 0.023 (equation (14)) multiplied by 
60“-*, where 60 is the length-diameter ratio of the tube used 
in obtaining the data upon which equation (14) is based. 

The data for modified Reynolds numbers below about 
10,000 show a trend with respect to length-diameter ratio 
and cannot be represented by a single line. 



(a) Unoorrected lor eOect of lengtlxllametcr ratio. 

(b) Corrected for effect of lenith.dlametcr ratio. 

Fiqubi 9.— Effect of ratio of length to diameter using modincd method of correlating heat 
tratisfer coefficients. Inconel and platinum tubes; Ia1ct.alr temperature, S3S° R; physical 
propertlet of air evaluated at surface temperature. 

The data shown in figure 9 for the platinum tube (L/D, 40) 
were obtained for surface temperatures up to 3050® R and 
surface-to-air temperature ratios up to 3.5: it is of interest 
that, for this extended range of temperature, the data arc 
correlated with negligible scatter by use of the modified 
Reynolds number. 

Effect of air temperature.^ — The data presented thus far 
were obtained with an inlet-air tempfirature of about 535° R. 
A wide range of sm-face temperature and corresponding 
surface-to-air temperature ratio was investigated. In order 
to determine whether the correlation applied for a wider 


HEAT TRANSFER AND FRICTION COEFFICIENTS FOR AIR FLOWING IN SMOOTH TUBES 


351 


range of air temperature, tests were made in which, the inlet- 
air temperature, and hence average air temperature, was 
varied. The results are shown in figm’e 10. 

Inspection of figure 10 indicates that the ordinate param- 
eter decreases as the inlet-air temperature and the corre- 
sponding average air temperature increase. For example, 
in figure 10 (a), at a Reynolds number of approximately 
100,000 and a surface temperature of approximately 1800° R, 
an increase in inlet-air temperature from 535° to 1460° R 
(with a corresponding decrease in Tt/T^ from 1.9 to 1.2) 
results in a reduction in the ordinate parameter of about 24 
percent. The ordinate parameter would decrease further 
if the inlet temperature were increased until T,jTi was 
approximately 1.0. At constant values of TJTt, the 
ordinate parameter varies approximately as temperature 
to the —0.2 power. 

As much as one-fourth of the observed separation' of the 
data with temperature level may be caused by experimental 
error in measuring air temperature at the high temperature 
levels. Inasmuch as there is a considerable variation in the 
reported values of thermal conductivity at high tempera- 
ture, as is evident in figure 3, the effect may be partly or 
entirely the result of using improperly extrapolated values 
of thermal conductivity. 

The values of conductivity from references 7 and 8, which 
have been used thus far in correlating the present data, 
vary appro.ximately as the 0.85 power of the absolute 
temperature. The values from reference 10 vary approxi- 
mately as the 0.73 power of the temperature. Inasmuch as 
the ordinate parameter of figure 10 varies inversely as 
conductiAuty to the 0.6 power, the conductivity would have 
to vary as temperature to the 0.85 power minus 0.2 divided 
by 0.6, or to the 0.5 power, in order to completely eliminate 
the effect of temperature level in figm^ 10 (a) and (b). 

Inspection of figures 5 and 6 indicates that the use of 
smaller values of thermal conductivity at high temperatures 
in correlation of the data would also result in a decrease 
of the effect of surface-to-air temperature ra'tio, so that 
correlation of the data according to eq^uation (15), in which 
the fluid physical properties and density are evaluated at the 
surface temperature, would result in an overcorrection 
of the data with increasing values of surface-to-air temper- 
ature ratio. A relation of the form of equation (15) can, 
however, be used for correlation of the data by evaluating 
the fluid physical properties and density at a reference 
temperature lower than the surface temperature. For 
example, if conductivity is assumed to vary as temperature 
to the 0.73 power, the effect of surface-to-air temperature 
ratio can be decreased by using a reference temperature 
To .75 midway between and T,. If conductivity is assumed 
to vary as the square root of temperature, the effect can 
be eliminated by using the film temperature Tf as the 
reference temperature. 

In order to illustrate the effect of using smaller values of 
thermal conductivity at high temperatures, the data of 



(a) Physical properties of air avaloatedat surface temperature; thermal condnctlvitr propo> 

t&DoI to Inconel tube; length-diameter ratio, 120. 

(b) Physical properties of air evaluated at snrCice temperature; thermal condnctlrlty' pcopor- 

tfonal to piat fniim tube; length-diameter ratio, 46. 

Figubs 10. — Effect of inlet-olr temperature modified method of correlatlaf beat-transfer 

coefficients. 


figure 10 (a) are replafted in figure 10 (c) using the conduc- 
tivity data of reference 10, and with the physical properties 
and density evaluated at the film temperature T^.n- The 
line representing equation (15), which for the lower reference 
temperature becomes 


034 Ap.o^Ma.75 Y ‘^ /A'\ 

«0.75 \ M0.7S / V *0.75 / 


is included for comparison. 
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(c) Physical properties of air evaluated at fllm temperature Ti.n: thermal conductivity 
proportional to Inconel tube. 



(dj Physical properties of air evaluated at fllm temperature 7V; thermal conductivity propor- 
tional to -^T', Inoonel tube. 

F iGURE 10.— Concluded. Effect of inlet-air temperature using modified method of correlating 
heat-transfer ooeifidents. 


Tko data of figure 10 (b), although uoL shown, would 
fall within tlie total separation of the data in figure 10 (c). 

Figure 10 (d) is included to illustrate the effect of using 
still lower values of thermal conductivity at high tempera- 
tures. The data are the same as in figures 10 (a) and 10 (c), 
but in figure 10 (d) the conductivity is assumed to vary as 
the square root of temperature and the physical properties 
and density are evaluated at the film temperature Tf. 
The line representing equation (16), which for the reference 
film temperature Tj becomes 

is included for comparison. 

In figure 10 (d) it is indicated that both the effects of air 
temperature and surface-to-air tompeiature ratio are olimi- 
natfed by use of the assumed conductivity-temperature re- 
lation and by evaluation of the physical properties and 
density at the film temperature. The data can he repre- 
sented \vith good accuracy by equation (15b). 

No justification of the arbitrarily assumed comluctivity- 
temperature relation shown in figure 3 can bo given other 
than that it results in the best correlation of the data reported 
heroin. 

Effect of heat extraction. — In reference 12, it is predicted 
that with heat extraction from the air the Nussclt number 
would increase as the heat-flux density increases, or as the 
surface-to-air temperature ratio decreases. It was further 
found that the effect could bo eliminated b}' evaluating the 
physical properties and density' at a reference temperature 
midway between Ti and T,. 

Tests were conducted to determine whether or not the 
effect of T,!Tt on correlation of heat-transfer coefficients 
was the same for heat extraction as for heat additipn. These 
results were obtained with an Inconel tube having a length- 
diameter ratio of 60 and a bellmouth entranee for a range 
of T,[Tt from Q.46 to 0.82. The data are correlated in the 
conventional manner in figure 11 using thermal conduotlvi- 
ties from references 7 and 8 with the physical properties 
evaluated at the average temperature Included is the 
line representing equation (10) modified for variation in 
LID, which becomes 


For an inlet-air temperature of 1460° R, a surface tem- 
perature of 1850° R, and a Reynolds number of 85,000, 
the ordinate parameter is about 15 percent below the refer- 
ence line. Comparison of figures 10 (a) and 10 (c) indicates 
that the use of the lower values of conductivity from refer- 
ence 10 reduces the effect of temperature level by approxi- 
mately one-half, or for constant v^ues of TJT^, the or^nate 
Iiarameter varies as temperature to the —0.1 power. 


The effect of a decrease in T./Tj, from 0.82 to 0.46 on the 
correlation of the data is negligible, and the data can bo 
adequately represented by equation (16). 

The data in figure 11 do not necessarily show that no 
effect of TJTi exists for cooling inasmuch as the use of a 
smaller variation of conductivity with temperature would 
result in an increase in Nusselt number for a decrease in 
T,ITi and hence introduce an effect of r,/T» into the data. 
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Reynolds number, 

Ficcee 11. — Oonventlonal method of conelathig heat-transfer coefflcfenta with heat e itrae 

tibn ftom air. Inconel tube; length-diameter ratio, 60; beTlmoiith entrance. 

The effects of using smaller variations of conductivity are 
shoivn in. figures 12 (a) and 12 (b) with the data based on 
the conductivity data of reference 10 and on the assumption 
that conductivity varies as the square root of temperature, 
respectively. In figure 12 (a), the date are correlated ac- 
cording to equation (ISa) for heat addition. A considerable 
overcorrection for the effects of T,[Tj, is obtained, which 
indicates that the effect of TJTt for cooling is smaller than 
that obtained for heating. In figure 12 (b), however, where 
the data are correlated according to equation (15b) for heat 
addition, the separation of the data is very small, mdica,ting 
that for the assumption that conductivity varies as the 
square root of temperature the effect of TJTf, is the same 
for both the heating and cooling data, and that both the 
cooling and the heating data can be represented by equa- 
tion (15b). 

■ Summary plots of heat-transfer data. — ^Heat-transfer data 
have been shown for a range of Reynolds number, surface 
temperature, inlet-air temperature, entrance configuration, 
and length-diameter ratio, and for heating and cooling of 
the air. The results are summarized in figure 13, which 
shows on single plots date from previous figures, which 
include aE the foregoing variables. 

The date are corrdated in figure 13 (a) bj" the modified 
method using the thermal conductivity data of references 7 
and S with the physical properties and density evaluated at 
Tf The line representing equation (15) is included for 
comparison. As previously shown, aE the heating date ob- 
tained with an inlet-air temperature of 535° R regardless of 
length-diameter ratio, surface temperature, or entrance con- 
figuration can be represented with good accuracy by the 
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(b) Physical properties of afr evalnated at film temperature Tr, thermal conductivity pro- 
portional to Vr. 

Figuee 12.— Modified, method of correlating heat-feanafer coetSoients with heat extraction 

from air. 


reference line. The high inlet-air temperature data and the 
cooling data faE below the reference line, so that equation 
(15) does not apply with desirable accuracy for aE the 
conditions investigated. 

The date of figure 13 (a) are replotted in figure 13 (b) 
using the thermal conductivity date from reference 10, and 
with the physical properties and density evaluated at To.rs. 


354 


REPORT 1020— NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



Modified Reynolds number^ 


(a) Physical properties of air evaluated at sor^oe temperature T«; thermal conductivity 
proportional to 



Modified Reynolds number, 


fh) Physical properties of air evaluated at film temperature T».n; thermal otmductivlty 
proportional to T*-'*. 

Fiouas 13,— Summary plot of experimental data using modified method of correlating hcat- 
transfer ooelfteients. 

The line reprosonting equation (15a) is included. Trends 
with temperature level and with T./Tj are reduced compared 
with figure 13 (a). The high inlet-air temperature data and 
the cooling data again deviate farthest from the reference 
line. 



<c) Physical properties of air evaluated at film temperature Tr, thermal conductivity propor- 
tional to -{t. 

Fiouhe 13.— Concluded. Summary plot of experimental data using modified method of 
correlating heat-transfer cocfllelents. 

Figure 13 (c) shows the data of figures 13 (a) and 13 (b) 
replotted using the assumption that thermal conductivity 
varies as the square root of temperature, and with the, 
physical properties and density ovaluateitl at Tf. The lino 
representing equation (15b) is included. The conductivity- 
temperature relation assumed in figure 13 (c) results in the 
best correlation of all the present data. Trends with respeet 
to temperature level and with T,ITt for both heat addition 
and heat extraction are almost completely eluninatcd and 
the data can be represented by equation (15b). 

COHBELATION OF FEICTION COEFFICIENTS 

Friction coefficients without heat transfer. —Average half- 
friction coefficients //2 calculateil by equation (6) for flow' 
without heat transfer are plotted against Reynolds number 
GDjui, in figure 14 (a). Included for oompanson is the line 
representing the K&rmfin-Nikuradse relation betw'eeu friction 
coefficient and Reynolds number for incompressible turbulent 
flow’^, which is 



and the line representing the equation for laminar flow in 
circular tubes, which is 

2 GD 

Mi 


( 18 ) 
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Reasonably good agreement with the Ktan&n-Nikuradse 
line is obtained for aH five configurations in the range of 
Reynolds number corresponding to turbulent flow. It is of 
interest that the K§nn4n-Nikuradse relation, which was de- 
rived for incompressible flow, applies mill good accuracy to 
compressible flow, inasmuch as the present data cover a range 
of tube-exit Mach number up to 1.0. 

Figure 14 (a) illustrates the agreement of the data obtained 
in the present investigation without heat transfer with the 
results of other investigators, which are in general agreement 
with equations (17) and (18). 

riiction coefficients with heat addition. — ^Representative 
friction data, obtained with heat addition to the air,areshown 
in figure 14 (b). Data for no heat transfer {T,ITt, 1.0), and 
the Karmin-Nikuradse and laminar flow lines are included 
for comparison. These data are for an Inconel tube having a 
length-diameter ratio of 60 and a bellmouth entrance, and 
were obtained simultaneously with the heat-transfer data of 
figure 5. 

The data for low surface temperatures {T„ 605° and 675° R) 
agree fairly well with those for no heat addition. As the 
surface temperature is increased, however, //2 decreases at 
high Reimolds numbers and stays the same or increases at 



(a) No heat transTer; Inconel tubes; inlet-air temperatine, <!39° B. 

(b) Heat addition; Inconel tube; length-diameter ratio, 60; bellmoutb entrance; inlet-air 
temperature, S36° E. 

Figubs H.— Voriatiou of average half-lrfctian coefficient with. Beynolds number. 



(a) Viscosity and density evaluated at surthce temperature T,. 

FtorsE 15.— Varbtion of modified half-hdctiou coefllctoit with modified Beynolds number 
with beat addition. Ineonel tube; lengtb-dinmeter ratio, 60; bellmouth entrance; Inlet-olr 
temperatifte, 535“ E. 

low RejTLolds numbers, so that the data cannot be repre- 
sented by a single line. These results are typical of those 
obtained with other entrances and length-diameter ratios. 

Modified Motion coefficient. — ^Various methods of defining 
the Reynolds number and the dynamic pressure on which to 
base the average Metion coefficient were tried in an attempt 
to eliminate the trends shown in figure 14 (b). One such 
method was to base the dynamic pressure on a density evalu- 
ated at the surface temperature and a velocity evaluated at 
the average air temperature and to plot the resulting friction 
coefficient against the modified Reynolds number ptYiDjn,, 
The Metion coefficient defined in the foregoing way f, is 
given by equation (9). 

The results of this method of correlating the Metion data 
are illustrated in figure 15 (a), where /,/2 is plotted against 
p,VtD/iis. The data are the same as in figure 14 (b), and 
curves representing the Kfirmfin-Nikuradse and the laminar 
flow equations are again included. 

In general, the Metion coefficient is overcorrected and 
increases with an increase in T,[Tt at constant modified 
Reynolds number. Figure 15 (a) is typical of the results 
obtained with other entrances and length-diameter ratios. 

The friction factor was also defined on the basis of a 
density evaluated at the film temperature Tf (equation (8)). 
The resulting Metion coefficient was then plotted against the 
modified Reynolds number with the density and viscosity 
evaluated at the film temperature Tf. The data of figure 15 (a) 
are replotted in this manner in figure 15 (b). The trends 
with T,[Tj, are eliminated for Reynolds munbers above 
20,000 and the data can be represented wdth reasonable 
accuracy by the Kfirmin-Nikuradse relation using modified 
parameters. 


356 


REPORT 1020 — ^NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Effect of entrance configuration. — Friction data obtained 
simultaneously with the licat-transfer data of figure 8 using 
an Inconel tube (length-diameter ratio, 60) with long- 
approach and right-angle-edge entrances are shown in figure 
16 (a). Data for no heat transfer are included. At Reynolds 
numbei's above 20,000, the data of figure 16 (a) are in reason- 
able agreement with the K^rm&n line, which was seen in 
figure 15 (b) to fit the belhnouth-entrance data \uth good 
accuracy. The data for the right-angle-cdge entrance fall 
shghtl^' below the reference line. No pai’ticular significance 
can bo attached to this small deviation, however, and it is 
felt that the average friction coefficients for all tlnve entrances 
can bo approximated with sufficient accuracy at modified 
Reynolds numbers gi'cater than about 20,000 by the relation 



Cb) Viscosity and density evaluated at film temperature 7>. 

FiauRA 16.— Concluded. Variation ot modified half-fricUon ooeffident with modified Key* 
Holds number with heat addition. Inconel tube; length-dlametor ratio, 60; bellmouth 
entrance; Inlet-alr temperature, 636^ B. 

Effect of length-diameter ratio and inlet-air tempera- 
ture. — Friction data obtained with heat addition to the air 
for length-diameter ratios of 30 and 120 and for inlet-air 
temperatures from 635° to 1460° R are shown in figure 16 (b). 
The data of figure 16' (b), which are for a bellmouth 
entrance, fall above the K&rm&n line, whereas the bellmouth 
data in figure 15 (b), which correspond to a length-diameter 
ratio of 60, fall on the line. No particular significance can be 
attached to these differences, however, and it is felt that the 
data for tlie range of length-diameter ratio and inlet-air 
temperature considered can be represented with sufficient 
accuracy by equation (19). 



(a) ESect of entnnoe conflguratton with heat addition. InooncI tube; Icnsthaiameter 

ratio, 60; inkt.alr temporature, S3S° R. 

(b) Effect of Icngth^iiameter ratio and Inlet-air temperaturo with heat addition. Inconel 

tubes. 

(o) Effect of heat extraction. Inconel tube; Icngth-.illamctor ratio, 60. 

FiauBi 16.— Variation of modified half-friction oocfllclcnt with modified ReynoMs number. 

Effect of heat extraction. - -Friction coefficients obtained 
simultaneously with the data of figure 1 1 for heat extraction 
from the air are shown in figure 16 (c). Data for no beat 
transfer arc included. The coordinates and the reference 
line arc tlie same as in figures 15 (b) to 16 (b). The data can 
be represented with reasonable accuracy by the referenee line 
for turbulent flow. 
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SUIVIMARY OP RESULTS 

The results of this iavestigation of heat transfer and 
associated pressure drops conducted with air flowing through 
smooth tubes for a range of surface-to-air temperature ratio 
from 0.46 to 3.5, surface temperature from 535° to 3050° R, 
inlet-air temperature from 535° to 1500° R, Reynolds num- 
ber up to 500,000, exit Mach number up to 1.0, length- 
diameter latio from 30 to 120, and for three entrance con- 
figurations may be summarized as follows: 

1. Gonventional methods of correlating average heat- 
transfer coefficients, whereiu specific heat, viscosity, and 
thermal conductivity are evaluated at the average air 
temperature, resulted in a considerable and progressive 
separation of the data with increased surface-to-air temper- 
ature ratio for heat addition to the air. Evaluation of the 
physical properties at the film or surface temperature, instead 
of the bulk temperature, aggravated the separation. 

2. Satisfactory correlation of heating data obtained with 
inlet-air temperatures near 535° R was obtained for the range 
of surface temperature considered when the Reynolds number 
was modified by substituting the product of air density evalu- 
ated at the surface temperature and velocity evaluated at the 
average total air temperature for the conventional mass flow 
per unit cross-sectional area, and the properties of the air 
were evaluated at the average surface temperature. 

3. Correlation of the heating data obtained at inlet-air 
temperatm-es above 535° R by the modified method resulted 
in a decrease in Nusselt number for an increase in average 
air temperature, which was eliminated by the assumption 
that the thermal conductivity of air varies as the square root 
of temperature and by the evaluation of the physical prop- 
erties of air including density at a temperature midway 
between air and surface temperatures instead of at the 
surface temperature. 

4. With tile assumption that thermal conductivity varies 
as the square root of temperature, data obtained for cooling 
of the air indicated that the effect of surface-to-air temper- 
ature ratio was the same as for heat addition, so that both 
heating and cooling data can be represented by the same 
equation. 

5. The use of different entrances to the teat section (that is, 
a long-approach, a right-angle-edge, or a beUmouth entrance) 
had negligible effect on average heat-transfer coefficients, at 
least for length-diameter ratios of 60 or greater. 

6. Yariation in tube l eng th-diameter ratio had a relatively 
small effect on average heat-transfer coefficients, which was 
satisfactorily accounted for by including a parameter for 
length-diameter ratio in the correlation equation. 

7. Friction coefficients for no heat transfer, calculated 
from a dynamic pressure based on an average air density, 
were in good agreement with those obtained by other inves- 
tigators. The corresponding friction coefficients obtained 


with heat addition to the air showed a considerable effect of 
surface-to-air temperature ratio for the Reynolds numbers in 
the turbulent region. 

8. All the friction data could be represented with reason- 
able accuracy by the K&rmSn-Nikuradse relation for incom- 
pressible turbiilent flow by using a friction coefficient cal- 
culated from a dynamic pressure based on a density evalu- 
ated at the film temperature and the modified Reynolds 
number in which the density and viscosity were evaluated at 
a temperature midway between the air and surface tem- 
peratures. 

Lewis Flight Propulsion Laboratoht 
■ National Advisort Committee for Aeron.autics 
Cleveland, Ohio, December Sl^ 1950 
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